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The detailed operation of an electron spin entangler is studied, using density matrix equations. 
The device is made of a superconductor, two quantum dots and two normal leads. The 
treatment takes into account coherent tunneling in a non-perturbative way, and analyzes the 
various parasitic effects, in addition to the main process (crossed Andreev reflection) : those 
include singlet pairs passing through a single dot, or cotunneling between dots through the 
superconductor. The optimum operation of the device is characterized. 



Producing entangled electronpairs is a challenge for fundamental experiments (analogous 
to those performed with photonP), as well as controlling quantum information^ in solid state 
deviceM It has been pro pose d that a superconductor, connected to energy filters, can serve 
as a source of spin singlets™. Here we address the operation of the "S-DD entangler" made 
of a superconductor (S) connected in parallel to a double quantum dot (DD) where Coulomb 
blockade prevents two electrons to pass through a single dot. The dots are s mall enough so that 
a single electron state is involved in each dot. The "crossed Andreev" (CA)^^ process is thus 
favoured, where a spin singlet is emitted, shared by the DD. The feasibility of such a device 
crucially depends on the control of "parasitic" processes spoiling entanglement, mainly of two 
kinds : first, singlet pairs can pass through a single dot, either through a double occupation 
state (direct Andreev DA), or one by one (Fig. 1). Second, an electron can pass from one dot to 
the other by elastic cotunnelin{piHI (CT) through S. The first, but not the second, was considered 
in Ref. 5. In addition, all processes are mixed together, making a consistent treatment difficult. 
Such a study is indeed possible by using the density matrix equations, which generalize the usual 
master equations to the inclusion of coherent processes. Those correspond to both Andreev 
transitions or to cotunneling. They are made of two virtual transitions, with a quasiparticle 
created in S then destroyed. On the contrary, single electron transitions between dots and the 
leads L,R are incoherent. The complete quantum master (QM) equations for the subsystem 
made of S and the two dots can be derivecPfor instance following Ref. 10. Here the discussion 
is based on the analysis of the dot populations and averaged current flow. Further results are 
devoted to sho t no ise correlations and Bell inequalities^! Notice that the S-DD entangler was 
recently studieJP^ in se ries with a splitter detecting entanglement! On the other hand, QM 
equations were employe cP for a different principle of entangler using another dot instead of a 
superconductoP^. Also, QM equations were a t th e basis of the analysis of a device permitting 
teleportation of the electron spin in a dot array^. 

In the ideal operation of the S-DD entangler, the Coulomb blockade in each dot is strong 
enough so as to rule out double occupancy. Starting from an empty state 00, CA reflection allows 
transitions to the sin glet state, shared between the two dots, ll s with a rate 7,4 T (7,4 is the 
geometrical factor^^). For a resonant CA process, the dot energy levels satify e = E\ + Ei = 0. 
The two electrons are evacuated from the dots into the reservoirs (with chemical potentials 
Hl,r < E\,E2) and the transitions to states 01, 10 occur with rates Tj (i = L,R). 

If the Coulomb charging energy is not so strong, a coherent transition from 00 to a doubly 
occupied dot states 20 or 02 can occur via a direct Andreev (DA) process, which has a rate Tj 
which is larger than for the CA process (Fig. 1). Electrons can subsequently be detected into 
this reservoir, with rates P^. This conduction channel implies dot energies Ei + Ui associated 



with double occupancy. One may also start from an initial state 10 or 01. DA can then 
proceed through the empty dot, but the charging energy of states 21 or 12 equals E\ + E2 + t/». 
Detection in the reservoir can either lead to 20 (02), or to singlet and triplet states ll Sj t. Another 
parasitic channel involves two electrons of a Cooper pair tunneling one by one towards the same 
reservoir (Fig. 1). This involves a singly occupied virtual state which costs an energy A5, the 
superconducting gap. Contrary to the DA process, the dot is emptied before the quasiparticle in 
S is anihilated. This process happens with a rate = I^Tf / 'A|. Last but not least, cotunneling 
(CT) allows a coherent transfer of an electron from one dot to the other via S. It couples states 
01 and 10, but also 20 (02) and 11, 21 and 12. CT has a rate "fcT which is reduced by a 
geometrical factor. 
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Figure 1: Full operation of the entangler, including the three Andreev channels and cotunneling. Real states are 
squared while virtual states are dashed squared. States with three electrons are omitted for clarity, spin is not 
represented. The Si's are the one-electron tunneling matrix elements. The crossed Andreev (CA) process is in 
the middle, the direct Andreev ones (DA) second from top or bottom, and the one-by-one processes are the top 
or bottom ones. Cotunneling (CT) connects the states vertically. 

Let us exclude high energy states 22 with N = 4 electrons in the double dot (DD), keeping 
states 00, aO, Oct (with spin a), 11 (singlet and triplets), 2a and a2 (2 means a local singlet). 
Transport of electron pairs through the DD is highly correlated. With CA alone, pairs pass 
one after the other through the DD. Including the parasitic processes, there is a strong mixing 
of processes and no simple perturbative calculation is possible. Yet, starting from the full 
one-electron Hamiltonian, one can derive a complete set of QM equations for the populations 
Pa = Paa and the "coherences" p a p (diagonal and non-diagonal matrix elements of the density 
matrix pdd)- It takes the general forrrJ^ : 
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i(E a — Ep)p a p + i^2(p ai Tp 1 - /5 7 /?r a7 ) — *^2(Ta-y + T/3 7 ) 



(1) 
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where the fi's are the coherent rates and the T's the incoherent rates. One assumes in the 
derivation that virtual states with at most one quasiparticle in S are created. The obtained 
set of equation^ is valid up to second order in the matrix elements describing tunneling to the 
leads, and to any order in the coherent rates (which are of order two in the tunneling matrix 
elements between the superconductor and the dots). 



In Ref. 5, a T-matrix calculation was performed, calculating separately the ideal (CA) 
current, and the DA and one-by-one parasitic currents. Here the optimum operation of the 
device can be settled on a firmer basis, and a better understanding of the physics involved is 
obtained. First let us assume a symmetric device (Tl = Tr) and treat the processes separately, 
without cotunneling. Assuming j A T 2 2> e 2 , the CA current in each lead is 



I C L A « eT-, 



8 7 2 T 2 



8 7 ^T 2 + T 2 /4 

On the other hand, the DA current I DA and the one-by-one current I obo read 

16T 2 ^ KT 2 T 



T DA -p 



lt° = e- 



16T 2 + r 2 + ?7 2 ' * L ^A 2 +KT 2 
where K is a numerical constant. If T « U, As, one has 
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(4) 



(5) 



The general case can be treated, setting e = E\ + E<i and putting all processes together. As 
an illustration, an analytical formula can be given for the total current, up to first order in the 
parasitic processes 



II = ea [T L + T R + AT 2 - ^ - 2KAa T L ^ (l - ^) 
-8Aa r L £ (5 - i) - 2a r L ^? ( 5 - i + r2+8 2 ^ +£2 )] 

where A = r2+ ^|^l +e2 and ^ 1 = A + 1 + iVr L + Ti/r^. 



(6) 
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Figure 2: Populations of the DD states as a function of U/F for As = 9.5K, E\ — —E2 = 0.5K, F — 
T = 0.1K, 7a, 7c = 0.2. States b, c, d, e, f, g, h, k respectively correspond to charge states 

(ll s ),(10),(01),(20),(02),(21),(12) and (11 ( ). 

The optimization of the entangler requires that most pairs emitted by the superconductor 
leave separated into L and R. To fight DA and one-by-one processes, one may compare the 
currents, taken individually for each process, or alternatively adopt a dynamical argument : 
starting from state \a) = (00) at time t = 0, the probability pn of singlet state \b) oscillates 
slowly (with frequency 7^4 T), but with a large amplitude. On the contrary, the probability of 
state |e) = (20) oscillates more rapidly (with frequency T), but with a small amplitude. The 
competition between the two processes crucially depends on the decay rate T. If it is small, CA 
is favoured, but if it is too large, DA process wins, state \b) has no time to form. A detailed 



analysis gives the criterion J7, A<; >> max\T,T /7a> £ /7a]- A similar analysis can be made for 
the effect of cotunneling : once in the state (10) or (01), decay in L,R must be faster than the 
cotunneling frequency 7c T, leading to the criterion ^qT « max[\Ei — E^T]. This can be 
confirmed by a numerical solution for the probabilities of various states (Fig. 2). 

In summary, density matrix (quantum master) equations can be derived from a microscopic 
Hamiltonian for a realistic entangler, and allow to integrate all processes in a coherent and 
non-perturbative way. A range of parameters for optimum operation is ^aT,^cT « Yl^r « 
U,As- More info rma tion can be obtained by further analysis of the current fluctuations (shot 
noise correlations^**^. 
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